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The study presents an analytic solution of the Balitsky-Kovchegov (BK) equation in a particular 
kinematics. The solution is written in the momentum space and based on the eigenfunctions of the 
truncated Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation in the gauge adjoint representation, 
which was used for calculation of the Regge (Mandelstam) cut contribution to the planar helicity 
amplitudes. We introduce an eigenfunction of the singlet BFKL equation constructed of the adjoint 
eigenfunction multiplied by a factor, which restores the dual conformal symmetry present in the 
adjoint and broken in the singlet BFKL equations. The proposed analytic BK solution correctly 
reproduces the initial condition and the high energy asymptotics of the scattering amplitude. 

PACS numbers: 12.38.Bx, 12.40.Nn, 11.55.Jy, 24.85.-|-p 


I. INTRODUCTION 

The power growth of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [l| amplitude contradicts the unitarity condition 
and should be modified. In the pioneering paper of Gribov, Levin and Ryskin (GLR) 0 the authors introduced a 
gluon saturation, meaning that the gluon density grows at large energy resulting in an overlap of the gluon states which 
tarns growth of the amplitude with increase of energy. This mechanism introduces a new scale commonly referred to 
as a saturation scale Qs- The GLR equation was first to account for non-linear effects for describing the physics of 
the Deep Inelastic Scattering (DIS). The GLR equation has a linear term representing BFKL type ladder diagrams, 
but also it includes a non-linear term, which stands for so-called fan diagrams describing the pomeron splitting into 
two pomerons. The GLR equation was formulated in double logarithmic kinematics, where one resums powers of 
logarithm of energy and typical transverse momentum. This kinematics is more restricted with respect to the BFKL 
multi-Regge kinematics, where the energy of the colliding particles is much larger than any transverse momentum. A 
decade after formulation of the GLR equation the unitarization corrections in the Regge kinematics were taken into 
account in the further generalization of the BFKL equation in the Wilson line formalism by Balitsky Q and soon after 
that independently by Kovchegov in the color dipole approach Q . The Balitsky chain reduced to the equation 
derived by Kovchegov in the limit of infinite number of colors. The resulting Balitsky-Kovchegov (BK) equation has 
a relatively simple form in the transverse coordinate space where the integral Kernel of the linear and the non-linear 
terms are happened to be the same and have a simple interpretation of splitting of the color dipole into two dipoles. 
Despite numerous attempts during almost two decades the full analytic solution to the BK equation is still to be 
found. The linear part of the BK equation corresponds to the BFKL evolution, which was completely solved exploiting 
conformal symmetry of the BFKL Kernel in the space of transverse coordinates. However, the direct substitution of 
the BFKL eigenfunctions in the coordinate space to the BK equation does not lead to significant simplifications in 
solving the BK equation. It was recently proposed that BFKL equation in color adjoint representation can be solved 
for non-zero transverse momentum using the eigenfunction in the momentum space. This was done for calculating 
the Regge (Mandelstam) cut contribution to the planar helicity amplitudes in Regge kinematics. A resulting series 
of publications confirmed this calculations to be correct by performing analytic continuation followed by Regge limit 
of the so-called remainder function of the Bern-Dixon-Smirnov (BDS) scattering amplitude. The analytic solution 
of the color adjoint BFKL was also useful for calculating Regge cut corrections to the planar amplitudes of the next- 
to-maximally helicity violating (NMHV) configuration of the scattered particles Q. The Moebius representation of 
the adjoint BFKL equation was calculated in ref. Q. 

In the present study we build a general form of the analytic solution to the BK equation (see eq. (flUl) ) based on the 
eigenfunctions of the color adjoint BFKL equation. The difference between the eigenfunctions of the adjoint BFKL in 
eq. (in]) and the eigenfunctions of the singlet BFKL in eq. dZD is mainly in the proper normalization condition leading 
to a factor, which is related to a shift the argument of the BFKL eigenvalue by 1/2 with respect the eigenvalue of 
adjoint BFKL. This factor is needed to restore the conformal symmetry in the momentum space of dual coordinates 
Pi = Zi — Zi+i, i = 1,2 present in the adjoint BFKL and broken in the singlet case. The dual conformal symmetry 
is believed to be a sign of integrability in planar amplitudes and is supposed to be present in the integrable singlet 
BFKL equation. 

Next, assuming a particular dependence of the impact factor on the transferred momentum in eq. dSOj) we find a 
closed analytic solution of the Balitsky-Kovchegov equation in this particular kinematics. Finally, we check that our 
solution correctly reproduces the initial condition and the high energy asymptotics of the scattering amplitude. Some 


2 


details of the calculations are presented in the Appendix. 


II. SOLUTION OF THE BK EQUATION 


The Balitsky-Kovchegov (BK) ii equation describes the energy evolution of the imaginary part of the scattering 
amplitude N (x, y) as follows 

(9y A(x, y)=^ y y) ■ y) ■ ^) n { z , y) ^. (i) 

The BK equation is written in terms of ‘t Hooft coupling 
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and the rapidity 
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where s is the square of the total energy of the scattering particles and sq is some energy scale. Numerous attempts 
of solving BK equation in coordinate space suggest us lookin g fo r analytic solution in other representations. It is 
natural to start with momentum space using a notation of ref. [10| for the corresponding Fourier transform 
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and write the Balitsky-Kovchegov equation in the momentum space as follows 
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Before considering the non-linear evolution we go back to the linear BFKL equation, which can be solved in two 
different ways; in the coordinate and in the momentum spaces. Both solutions are well known and usually addressed 
to as q = 0 solution (momentum space) and q ^0 solution (coordinate space). Those two were thoroughly discussed 
in numerous studies and we do not focus on them for the sake of brevity. For more details on these two solutions 
the reader is referred to review texts [HI- This is true for the BFKL equation projected on the color singlet state of 
two interacting reggeized gluons. However, if one considers color adjoint BFKL with some infra-red divergent parts 
removed, there is another q 0 solution in the momentum space expressed through the adjoint BFKL eigenfunctions^ 


f 1^,71 (fc; q) 
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which were used in a serious of publications [l^ for calculations of the Regge (Mandelstam) cut contribution to 
the helicity amplitudes in maximally supersymmetric theory. Its solution was written in terms of f,y,nik,q) and 
a corresponding eigenvalue of the adjoint BFKL equation of shifted argument with respect to the singlet BFKL 
eigenvalue. The corresponding adjoint BFKL Green function was investigated in ref. [l^. 

In the present study we consider a non-linear generalization of the linear BFKL equation, namely the Balitsky- 
Kovchegov (BK) evolution equation. We introduce an ansatz for the BK solution in a particular kinematic regime, 
which can be written in terms of functions similar to fu,n{k,q)- We slightly redefine the ’’adjoint” fu^n{k,q) to 
account for a proper normalization of the singlet BFKL as well as a minus sign, which happens to be important for 
the non-linear term of the BK equation as follows 
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^ Here we use complex coordinates defined hy k = kx iky 
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The function in eq. © solves the leading order singlet BFKL equation as we show below. The major difference 
between the adjoint and singlet eigenfunctions is the factor |^||^_^ |, which restores the dual conformal symmetry 
of the singlet BFKL equation. This analysis is beyond the scope of the present study and will be published by us 
elsewhere. 

The corresponding orthogonality condition is given by 


d^kF:Jk,q)F,. 


'{k,q) 


S(iy 




and the completeness condition reads 



di' K,n(k',q) F^,n(k,q) 
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Then we can write the analytic solution of the Balitsky-Kovchegov equation in the momentum space as follows 


A^(k,q)= E 


dz/ Cu^-yiiF') Fj^^^i^kj q^^ 


( 10 ) 


where the coefficient function Ci,^n{X) is found below plugging eq. (ITUl) in the BK equation. Here we make a basic 
assumption that the coefficient function C^^niX) does not depend on k or q. We discuss this approximation in more 
details below, where we consider the initial condition and the high energy limit of the obtained BK solution. 

It is useful to introduce new complex variables 


k 

k — q 


= w, 



In terms of the coordinates w and w* the function in eq. 0 reads 
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and the completeness relation in eq. is given by 
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Consider the first linear term of the BK equation in the momentum space eq. © 
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In a similar way the second linear term of the BK equation in eq. m can be written as 
1 ds 
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and finally the third linear term reads 
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The sum of these three linear terms gives an expression, which determines the BFKL eigenvalue 
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through a function 
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which is expressed in terms of the digamma function defined as a logarithmic derivative of Euler gamma function 
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In deriving eq. (1^(71) we used the identity 
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For more details on the calculation of ^) fke reader is referred to review texts 0. This way we have shown that 
Fv,n{kiq) in eq. ([7]) is indeed the BFKL eigenfunction with the eigenvalue u}{v,n) given in commonly used notation 
without any possible shift of the argument of digamma functions. 

Now we are in position of plugging F^ ^ik, q) into the BK equation in eq. ([5]). Namely, consider the non-linear term 
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where is calculated as follows (here we define k/k* = e*^"^) 
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Using the last expression we write 
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Note that the coefficient of the non-linear term typically associated with the triple Pomeron vertex, in this notation, is 
a number independent of anomalous dimension and conformal spin. Therefore the (z/, n) dependence of the BK solution 
enters only through eigenfunctions F^^n{k,q) and the BFKL eigenvalue a;(z/, n) in eq. d?]) and eq. (fT^ respectively. 

Using the expression for the non-linear term eq. (HH) in the BK equation we get an equation for the coefficient 
function C^^niX) as follows 


Its solution reads 


dXnX) 

dY 


aXniY) + bXn{yf- 


X.nX) 


Qa.n(O) 

a + b C'jy,„(0) (1 - e“^) ’ 


(27) 
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where a = oj{v,n) is the BFKL eigenvalue in eq. (HU) and b = —as/{\/2) is the coefficient of the non-linear term in 
eq. (1261) ^. The initial condition of the coefficient function C^^n(Y) is set by the following expression 


X,n(0) = \q? / dX' FXik',q) Xk',q) = \q\^ 


|w'- 11 




(29) 


The function $(fc,q) is some impact factor determining an initial condition of the non-linear evolution. The self- 
consistency of the approach suggests a choice of the impact factor describing multiple rescatterings, e.g. an impact 
factor of the Glauber-Mueller type (see also ref. [l^). Our basic assumption is that the coefficient function C'y_„(0) 
does not depend on the transferred momentum q. This means that we choose the impact factor in such a way that 


^{w',q) 




(30) 


in accordance with eq. (ESI) and eq. (USD, where ^(w) is some function of w and independent of q. This approximation 
seems to be reasonable for small values of the transferred momentum. 

The evolution equation eq. (1271) and its solution in eq. (1281) are similar to that of the phenomenological summation 
of the Pomeron Fan diagrams in the Reggen Field Theory for zero transverse momentum [l^ and can be used for 
summation of pomeron loops [ 13 . We recall that we wrote a general form of the BK solution in eq. OT as 

^ nOO 

Ar(k,q)= ^ / di^Xn{Y)Xn{k,q), (31) 

n = — rv-i ^ 


and it follows directly from this definition and the completeness condition in eq. (E]) that at K = 0 for our choice of 
Cv,nX) we have 


Ar(k,q)|^=0 = ‘J>(k,q). (32) 

A general form of the solution A/’(k, q) in eq. (1^ together with the expression for the coefficient function Cv,nX) 
in eq. ([281) and its initial condition (7i/,„(0) in eq. (l29l) present the main result of the this study. 

It is easy to see from eq. cni) and eq. ((^511 that the solution of the BFKL equation in the color singlet state is 
obtained by setting 6 = 0 in the coefficient function Cv,nX)- Note that in the BFKL equation, which corresponds to 
the linear part of eq. ([5]) the transferred momentum q is not mixed with the integration variable of the Kernel and 
thus the condition that C^^nX) does not depend on q can be relaxed reproducing the full analytic solution of the 
BFKL equation for a non-zero transferred momentum. 

We have shown above that our choice of functions F^n{k,q) is consistent with a solution of the BFKL equation, 
and the next step is to check that the full solution written in terms of these functions have a reasonable behaviour at 
high energy. 

In the limit K —>■ oo eq. (l28)) reduced to 


^ = V2 xXn) 


(33) 


^ Note that this coefRcient depends on the normalization of the eigenfunctions. 
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and the amplitude A/'(k, q) defined by eq. (1101) in this limit does not depend of the initial condition resulting in 

/•OO 2 

A/'(k,q)|i-_>oo - X! J dv V2 F^^ri{k,q) =—. (34) 


n——oo 


Details of this calculation are presented in the Appendix. This function peaks at small q and thus can be reasonably 
approximated by Dirac <5-function. 


Ar(k, q)\Y^^ 2<5(|gn<5(</,) = ^d(|g|)d(<^) = ^^(q). (35) 

Note that we have no k dependence at 1 " — > oo and the asymptotic amplitude diverges for q = 0. This amplitude 
represents so called fixed point at which the derivative dyJ^ik, q) vanishes and the amplitude goes to its maximal 
value. On the other hand one can search fixed points of the BK equation in the coordinate space in eq. o and it is 
widely accepted that a natural fixed point of eq. o is N(x.,y) = 1, which is also a fixed point in phenomenological 
models. Then for our definition of the Fourier transform in eq. we have 




(36) 


We have a certain freedom in choosing the cut-off c, which we fix by imposing a condition that the expressions in 
eq. (1551) and eq. (1551) are to be identical. This condition imposes a definite choice of the cut-off such that 


c2 = 


4g2V-(l)-2 

F ■ 


(37) 


Note that this choice of the cut-off can be related to two basic values that govern the saturation physics. Namely, the 
coefficient that determines the rapidity dependence of the saturation scale^ 


Q,(y)~g,(0)e“«2>', C7~4.88 

and the critical exponent 7 c = 0.6275 in 

Af(x,y) oc |x-yp'^'= ln|x-y|, for |x - y|0. 

Those two can be expressed at one percent accuracy through our choice of the cut-off as follows 

2 


V c^k'^ 


and 


In 


\ J 


-4.84 


- 0.634 


establishing a relation 


InC 


- 7c 


(38) 

(39) 

(40) 

(41) 

(42) 


The authors are not aware about any previous study discussing a similar relation. 


^ For most recent review on saturation physics the reader is refrerred to the book of Yu. V. Kovchegov and E. M. Levin and other review 
texts [Tl|l . 











7 


III. DISCUSSIONS 

We discuss a new form of the singlet BFKL eigenfunction in eq. ([7]) and using it build an analytic solution of the 
BK equation in eq. for a particular kinematics, where the impact factor depends on the transferred momentum 
through eq. dsnu- We check that the obtained solution satisfies the initial condition, it has a proper high energy 
behavior compatible with known values of the saturation scale and its form resembles a solution to phenomenological 
fan diagrams resummation. 

The Balitsky-Kovchegov equation was extensively studied during past two decades both numerically and analytically 
in various approximations. The comparison of the solution presented in this paper to the approximate solutions 
available in the literature is not a simple task and definitely requires a separate analysis which will be published 
by us elsewhere. The most important cross check for us is the compatibility with the linear evolution of the BFKL 
equation and the proper high energy behaviour as shown above. We expect some peculiar technicalities in going back 
to the coordinate space we have started with in eq. o, related to the fact that we are dealing with individually 
divergent quantities, which nevertheless result into a finite and well defined final answer. The proposed analytic 
solution of the BK equation may have many phenomenological applications [l^, e.g. effective high energy description 
of the proton-nucleus scattering as well as nucleus-nucleus scattering calculated through pomeron loops built of two 
symmetric (projectile-target) BK solutions. 

The singlet BFKL equation is reach of symmetries [l9l - l^ and being integrable is supposed to enjoy also the dual 
conformal symmetry as mentioned above. A similar though a more general dual superconformal conformal symmetry 
was identified for planar scattering amplitudes in maximally supersymmetric theory (see recent results in ref. |24l | 
and references wherein). The solution of the Balitsky-Kovchegov equation presented here is written in the space of 
transverse momenta and one can expect some manifistations of the dual superconformal symmetry appearing also in 
the case of the non-linear BK evolution. We leave these issues for our further publications. 
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Appendix A: Calculation of A7(k,q)|i 


We calculate the high energy asymptotics of the solution in eq. m 

OO «00 

A/'(k,q)|y_,.oo ^ / dv V2x(v,n) F^^ri(k,q), 

where the function x(i^, n) is given in eq. (I20[) by 
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Consider the following expression 
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where we use w/w* = We do the v integration first for |w| > 1 closing integration contour in the upper complex 
semiplane and the contribution comes only from one of the digamma functions in xi^-: Namely, 
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In a similar way we close the integration contour for v in the lower complex semiplane for |r(;| < 1 and get the same 
result. Using the definition of Fi,^nik,q) in eq (ITil) we finally obtain 


A/'(k,q)|y^oo 



dv V2 x(i 2 , n) F^^n{k,q) 
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